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Abstract. We study structural stability of smoothness of the maximal solution to 
the geometric eikonal equation on {M.'^,G), d > 2. This is within the framework of 
order zero metrics G. For a subclass we show existence, stability as well as precise 
asymptotics for derivatives of the solution. These results are applicable for examples 
from Schrodinger operator theory. 



1. Introduction and results 

In this paper we investigate the existence of a smooth global solution to the geometric 
eikonal equation on the Riemannian manifold (M'^, G), d>2, for a class of metrics G. 
We are interested in the so-called maximal solution S{x) constructed as the geodesic 
distance from x to a given fixed point xo (taken to be xq = 0). It is well-known 
that for some G"s this function 5* is smooth while for others this is not the case, in 
fact even S E might be false. (At this point the reader may consult [Lll ICC] 
for studies of "generalized solutions" to related Hamilton- Jacobi equations and further 
references.) Whence it is interesting to investigate the stability of smoothness of S under 
perturbation of the metric. Other issues we will study are bounds and the asymptotic 
behaviour of derivatives of a smooth 5* at infinity. This is a general mathematical 
problem motivated by specific applications in scattering theory, see |ACHt IBat HSj ISk] . 
Most likely it is relevant for other specific problems too (possibly from control theory, 
geometric optics, etc) although this will not be examined in this paper. 

To motivate our setup explained in further details below let us imagine a more general 
situation: Let us consider a complete simply connected (i- dimensional manifold {M,g), 
d >2, and a point Xq € M for which the exponential map $ = exp^^^{l-) : TM^^ — )■ M'^ 
is a diffeomorphism. Then the puUback G = is a metric on TM^g that can be 
compared with the canonical one, g^co- Introducing orthonormal coordinates we may 
identify TM^^ = and gxo{y,y) = (the usual Euclidean metric). Upon doing this 
identification the Gauss lemma (see |Cht Theorem 1.8]) implies that G{x) (considered 
as a matrix) has the form 

G{x) = P + P^G{x)P^, (1.1) 

where P denotes, in the Dirac notation, the orthogonal projection P = P{x) = \x){x\ 
parallel to 5: = x/\x\ and P± = P±{x) = I — P the orthogonal projection onto {x}-*-. 
Note that in this picture S{x) = \x\ and G(0) = /. 

Due to the above change of framework we can in principle reduce the study of a 
general stability problem to the one indicated above, i.e. for (R'^, G) only, in fact 
for the unperturbed metric being of the form (11. ip . In this paper we shall consider 
families {G{-)} that are of order zero, see (11. 2p and (11.30 below for precise definition. 
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Notice that this class is naturally equipped with norms giving precise meaning to the 
notion of "perturbation". Our main result asserts that, under conditions, indeed upon 
perturbing a metric of the form (11 .ip one obtains metrics, say denoted by Ge{-), for 
which the geodesic distance to the origin, S'e(-), is smooth (more precisely of class 
depending on conditions) and solves the eikonal equation 

VS,G;^VS, = 1 for a; ^ 0. 

Moreover introducing 

s^{x) = Se{x)/\x\ — 1 for a; 7^ 

we have bounds 

sup |5"s,(x)| = o(e°) for |a| < 2. 

a;eM''\{0} 

Depending on conditions there are somewhat similar bounds for higher order derivatives. 

Finally we use the change of frame in terms of the exponential mapping for the 
unperturbed metric, as explained above, to solve eikonal equations of the form 

\VS,{x)f = 2(A - V,{x)) for xeR'^\ {0}, 

where K is given by perturbing a negative radial function (potential) obeying certain 
properties, see Section [7l Here the parameter A > plays in applications (Schrodinger 
operator theory) the role of energy. The conventional way af constructing solutions is 
by a fixed point method, cf . for example \DS\ \B.61\ lis] . However to our knowledge this 
is not doable for our examples. 

This work is inspired by |Ba] which has similar results for perturbation of the special 
case G = I. As in [Baj we shall use the geometric/variational approach to define 
the maximal solution to the eikonal equation, see (II. 5p . Many of our arguments are 
however very different from those of |Ba] . This came out of necessity to treat the present 
generality, see Remarks II. 5l l3|). 1 1.71 and [5. II for comments on this issue. 

1.1. Conditions and main results. Let i5rf(M), ci > 2, be the space oidxd symmetric 
matrices with components in M and, for / > 2, let B^{W^) be the space of functions 
G:W^ ^ cSd(M) such that 

lldli = sup{(a;)l"l |9"c/ij(a;)| : x G MMa| < /, i, j = 1, . . . , 4 < cx), (1.2) 

where G{x) = (gij{x)) and (x) = (1 + Ixp)"*^/^. The space B^(W^) endowed with the 
norm defined in f ll.2p is a Banach space. Let M = MCMJ^) be the set of G G i3'(]R'^) for 
which there are positive constants a and b satisfying 

a\y\^<yGix)y<b\y\\ x,y eR''. (1.3) 

The set Ai is open in i5'(]R'^), and its elements will be referred to as metrics of order 
zero. We denote by Ti the Sobolev space {Hq{0, 1)Y with the norm 

\\hr = {h,h)= f\h{s)\'ds. 

Jo 
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For G G we consider the energy functional : M'^ x — ?■ M given, for {x, k) G M'^ x 7/ 
and y{s) = sx + k{s), by 

E{x,k)= [ y{s)G{y{s))y{s)ds. (1.4) 

Define for any G G a non-negative function S by 

S\x) = M{E{x, K):KeH};xe (1.5) 

where E is as in fll.4p . 
Let 

S^ = {ye {H\0, l)y : y{s) = sx + k{s), i^en}; xe M'^. (1.6) 

In agreement with usual convention [Mi] , we say that 7 G £x is a geodesic of G 
emanating from with value x at time one if di^E{x^ k) = as an element of the dual 
space %' of %. Geodesies can emanate from other points in W^. The general definition 
of a geodesic of G ( |Ch| IMij ) may be taken to be an orbit s — )■ 7(5) G M*^ solving a 
certain second order differential equation, see fl2.12p . These solutions can be extended 
to globally defined solutions, cf. the Hopf Rinow theorem |Cht Theorem 1.10]. By the 
Riesz lemma we can identify "H' and "H, as well as the set of bounded quadratic forms 
on "H by the set of bounded self-adjoint operators on "H. 

We shall consider two sets of metrics of order zero. The first one is the following. 

Definition 1.1. Let U be the subclass of order zero metrics G that satisfy: 

1) For every x G M'^ the metric G has a unique geodesic 7a;(s) = sx + Hx{.s) in E^. 

2) There exists c > independent of x G M*^ such that ii n = is given as in [1]) 

{dlE{x,K)h,h)>c\\h\\'^, heU. (1.7) 

Remarks. 1) The Hessian appearing in (II. 7p is given by 

{dlE{x,K)h,h2) 

= [ {2hiGh2 + 2yVG ■hih2 + 2yVG-h2hi+y{V^G-hi,h2)y)ds- (1.8) 
Jo 

here VG ■ h stands for the matrix {Wgij ■ h). Let us for completeness of presentation 
remark that for all G G satisfying [1]) there exists G > independent of x G M*^ 
such that 

\{dlE{x,n)h,h)\<C\\h\\'', hen. 

(This follows readily from ([ISD, and fIXTTiD .) 
2) For all G G W the corresponding exponential map $ = expo(l-) : TR^ ^ M'^ is a 
diffeomorphism, cf. [Cht Theorem 2.16] and [Mit Theorem 14.1]. 

Proposition 1.2. Let G ElA. The non-negative function S defined by (11.51) is of class 
G' on M'^ \ {0} and satisfies the eikonal equation 

VSG-^VS = 1 for Xy^O. (1.9) 

Furthermore, there exists G > such that 

sup(x)°>'"(l°l-^'l°l/2) |9"5(x)| < G for all \a\ < I. (1.10) 

|x|>l 
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Remark. One might suspect that (ll.lOp can be replaced by the stronger bounds 

sup(a;)l"l-^ < C for all |a| < I. (1.11) 

|x|>l 

In general this is an open problem. On the other hand f ll.lOp appears natural if the 
class of metrics is enlarged by replacing (11. 2p by 

= sup{(s)"^''^(l"l'i+l°l/2) Id^gijix)] : x G W^,\a\ < = 1, . . . , rf} < oo. 

In this situation indeed the analogue version of Proposition 11.21 holds true (showned by 
the same proof). 

Our second set of metrics of order zero is given as follows. 

Definition 1.3. Let O be the subset of order zero metrics G obeying: 

1) For X ^ 

G{x) = P + P^G{x)P^, (1.12) 

where P denotes, in the Dirac notation, the orthogonal projection P = P{uj) = 
{u\ parallel to u = x = x/\x\ and P± = -Pl(w) = I — P the orthogonal projection 
onto {uj}^. 

2) There exists c > such that for x ^ 

P±{G{x) + 2~^x ■ VG{x))P^ > cP^G{x)P^. (1.13) 

Note that since G is continuous at x = Definition ll.3| [T]) implies that G{0) = I. 
Moreover it follows from (I1.13P that c g]0, 1]. The simplest example of a metric of order 
zero satisfying (I1.12p and (11.130 is G = /, the d x d identity matrix. In Section [7] we 
provide other examples of metrics of order zero satisfying (I1.12p and (I1.13p . We may 
refer to (I1.13P as a convexity property since, given the orthogonal decomposition (I1.12p 
the estimate is equivalent to the geometric Hessian bound 

V'^Six^ > 2cg{x). (1.14) 

Here we use the conventional metric notation g rather than the matrix notation G 
and (with (I1.12p ) S{x) = \x\. Written in this way the condition (I1.13P clearly becomes 
geometrically invariant which a priori is a desirable property. However in computations 
we shall only use (I1.13p . 

In terms of the non- negative function S given by (II. 5p (for any G E M.) let 

s{x) = S{x)/\x\ — 1 for a; 7^ 0. 

Our first main result is: 

Theorem 1.4. LetU,0 Ai be given by Definitions \1.1\ and \1.3[ There exists a 
neighbourhood O Ai of O such that: 

i) O ^U; that is the set O is a subset of the interior ofU. 
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ii) Let G E O be given. Then there exist eo,C > such that for all G E M. with 
\\G — G\\i < Co not only G E O but also 

sup < CIIG-GII,, (1.15a) 
xeKd\{o} 

sup \x\\d''s{x)\ <G\\G -G\\f^ for\a\ = l, (1.15b) 

sup \x\^\d''s{x)\<G\\G-G\\y^ for\a\ = 2. (1.15c) 

Remarks 1.5. 1) One might presume that L{ is open in Ai. However if true at all this 
is a hard problem. Similarly (seemingly a softer problem) one might presume that 
if we drop the condition (11.131) of the Definition 11.31 and call this bigger class Oi, 
then Oi is a subset of the interior of U. (Note that Oi U due to Lemma 4.1.). 
Even this problem appears to be difficult. In our approach we use (11.131) crucially to 
obtain good control of perturbed geodesies uniformly in x. Given the lower bound 
in (II. 3p the condition (I1.13P is a somewhat weak assumption. 

2) The two constants eo, C > can be taken as locally bounded functions of (HGUi, a, c) G 
where the entries \\G\\i, a, and c are defined by (11.21) . (11.31) and (I1.13P for the metric 
G, respectively. Obviously this statement for / > 3 follows from the assertion for 
1 = 2. 

3) For perturbations of the Euclidean metric G = I one can replace the powers to the 
right in (I1.15ap - (ll.l5cp by the more natural factor ||G — G\\i, see |Baj . However the 
techniques of |Baj are not applicable in our case, see Remark 15.1 1 for some elaboration 
at this point. We do not know if this improvement is possible in our more general 
case. 

The estimate (11.15bp is a consequence, by interpolation, of the bounds (I1.15ap and 
(I1.15cp . cf. |Ho21 proof of Lemma 7.7.2]. We shall use (I1.15ap and the following weaker 
version of (11.15bp in the proof of (I1.15cp : 

sup |x| |9"s(x)| < CIIG-Gli;/^ for |a| = 1. (1.16) 

xeiR'*\{o} 

Our second main result supplements Theorem ll.4| [il|). It reads 

Theorem 1.6. Suppose I > 3. Let G E O and r > be given. Then there exist 
eo, C > such that for all G E M. with \\G — G\\i < and \a\ < I 

sup(x)l°l-i |9"^(x)| < C, (1.17a) 

|x|>r 

sup \d''S\x)\ < G. (1.17b) 

Remark 1.7. If / > 4 we can use Theorems 11.41 and 11.61 and interpolation to show that 
sup |9"s(x)| < Gr\\G - for 2 < |a| < / - 1. (1.18) 

|a;|>r 

A slightly improved bound for |a| = 3 was proved in [Ba] (in the setting of |Baj ) under 
the assumption that / = 3. We do not need estimates like (ll.lSp for |a| > 2 in our 
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applications [ISl ISk] . On the other hand Theorem 11.61 for / = 3 is indeed important in 
these apphcations. 

This paper is organized as follows: In Section |2] we study the minimization prob- 
lem (II. 5p . Using standard arguments we show the existence of a minimizer and some 
basic properties of any such minimizer. In Section |3] we show Proposition II. 2[ The 
proof is based on the implicit function theorem, and the somewhat lengthy scheme for 
proving the bounds (11.101) is used again in Section [6] to establish the improved bounds 
of Theorem 11.61 The proof of Theorem ll.4|[il) , given in Section HJ is based on an analysis 
yielding dynamical control of perturbed geodesies. We obtain sufficient control to be 
able to deduce the uniqueness of the energy minimizer from a result from global anal- 
ysis. In Section |5] we show Theorem ll.4| [ii|) by using various explicit computations in 
combination with results from Section HI The proof of Theorem II. 6[ given in Section 
[6l is based on functional analysis arguments for Hardy spaces taylored to the problem 
in hand. The first part of the proof is devoted entirely to setting this up abstractly. 
The second part is devoted to verification of conditions, and as indicated above, this 
involves a scheme from Section 121 In Section [7] we present examples from Schrodinger 
operator theory. 

2. The minimization problem 

In this section we study some basic properties for metrics G & Ai. Since S{x) = \x\ 
when G{x) = I then, from (11.31) and (II. 5p . we have 

a\x\^ < S\x) < b\x\\ (2.1) 

for all X eR^. 

Lemma 2.1. Let G E Ai and E as in (11.41) . Then for every x there exists k eH 

such that 

S\x) = E{x,k). (2.2) 
Moreover, if 'j is a geodesic of G emanating from with value x at time one then 

7(s)G(7(s))7(s) = /' 7(^)^(7^)7^^^, s E [0, 1]. (2.3) 
Jo 

In particular, if 'j is a minimizer of (ll.4p then for all s E [0, 1] we have 

7(s)G(7(s))7(s) = S'{x), (2.4) 

and rl^ajP < |7(s)P < -Isxl^ (2.5) 
b a 

Proof. Existence of a minimizer. To establish the existence of a minimizer it suffices 
to show that, for any fixed x E M.'^, the functional E{x, ■) is weakly lower semicontinuous 
on T-L, and that 

lim E{x, k) = oo. (2.6) 

IIkII^oo 



and 



-r\X\ 



<\7{s)\' <-\x\ 
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Note that (12. 6 p follows from the estimate, y{s) := sx + k{s) G £x, 

E{x,k)= y{s)G{y{s))y{s)ds > a \y{s)\'^ds. 



To prove the weak lower semi continuity of E in the second variable we let be 
a sequence in T-L that converges weakly to n E Ti and write y{s) = sx + k{s) and 
yn{s) = SX + Hn{s). Then 



E{x,Hn) - E{x,k) = / {yn-y)G{y){ijn-y)ds -2 j yG{y){ij - yn)ds 

Jo Jo 



yn{G{y) - G{yn))ynds. 







Since the first integral is non-negative, the second goes to zero as n — t- oo (because 
Vn — y converges weekly to zero) and the third integral goes to zero too (by the compact 
embedding of T-L into (C(0, 1))'^) we find that 



liminf E{x, k„) > E(x, k). 

n—^oo 



The conservation law. Since G is of class C' so is E{x, k) and for all h E H we 

have ^ 

{d^E{x,K),h) = I [yVG ■hy + 2hGy]ds. (2.7) 
Jo 

Moreover, if 7 is a geodesic of G then for all /i G "H 



f i'jVG ■ + 2hG'j)ds = 0. {2.i 
Jo 



Integrating by parts we obtain 



[ [7VG ■ hj - 2h{VG ■ 77 + G^)]dt = 0. 
Jo 



(2.9) 



Thus, in the space T-L' we have 

2efcG7 = 7VG ■ ekj - 2euVG ■ 77 (2. 10) 

for each k = 1, . . . ,d, where is the k^^ element of the standard basis of W^; ei = 
(1, 0, . . . , 0), 62 = (0, 1, . . . , 0), etc. This implies that 

^ = lG-^w, (2.11a) 

where the k^^ component Wk of the function w : [0, 1] — t- is given by 

Wfc = ■ efc7 - 2efcVG ■ 77. (2.11b) 

Letting e, denote the ordered basis ei, . . . , of M'^ we can write (12.1 lap and (I2.11bp 
more compactly as 

7 = K = 2-^G'"^7VG' ■ e.7 - G~^VG ■ 77. (2.12) 
Multiplying (I2.10p by % and adding over k yields 

jVG ■ 77 - 27VG ■ 77 - 27G7 = 0, 
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from which we obtain 

^(7(3)G(7(5))7(3)) = 0, .G[0,1], (2.13) 

and therefore (Q and (El) follow. 

The estimates. Suppose 7 is a minimizer of (II. 4p . Combining (II. 3p . (12. ip and (12. 4p 

we obtain the estimates for |7(s)| stated in (12. 5p . 

Furthermore, using the upper bound for |7(s)| we find that 

|7(s)| < ^ \i{r)\dT < ^J^\sx\, for all s G [0, 1]. 

Suppose now that for some si G (0,1] we have |7(si)p < (a/6)|sixp, then defining 
y{s) = (s/si)7(si) for < s < Si and y{s) = 7(5) for Si < s < 1 we have 

Jo Jsi 

Thus, using (11.31) . (12. ip . (12. 4p . and our assumption on si we find that 

E{y) < M^il^ + (1 _ s,)S^{x) < a\x\hi + (1 - Si)S\x) < S^{x), 
which is impossible by the choice of 7. Thus (12. 5p holds and the proof is complete. 

□ 

3. Proof of Proposition 11.21 
Clearly Proposition 11.21 follows from 
Proposition 3.1. Let G E U and write the unique geodesic for G with endpoint x as 

J^{s) = SX + K^.(s). 

i) The map M'^ 3 x — > /tx> G "H of class C^~^ . 
a) The non-negative function S in i\1.5\) obeys 

S\x)= [\,{s)GMs))Us)ds, (3.1) 

and it is a C' solution to the eikonal equation (11.90 . 
Hi) There are bounds 

< C'(x)i/2-"^'"(l"l-^/2>l/2) for all \a\<l- 1. (3.2) 

iv ) For any r > 

sup (a;)^i"(l°l-i' Id^'Six) I < Cr for all \a\ < 1. (3.3) 

|a::|>r 

v) In\m\) and\M) the constants C and Cr can be taken as locally bounded functions of 
{\\G\\i,a,c) G M.^ and {\\G\\i,a,c,r) G M^, respectively. Here the entries \\G\\i, a 
and c are defined by (11.21) . (II. 3p and (11.71) . respectively. 
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Proof. Re [i]). The statement follows from our assumption that G & U, the represen- 
tation ( \2.7\\ and the implicit function theorem. Note that indeed the classical implicit 
function theorem given for example in |De[ Theorem 15.1] (see possibly also Theo- 
rem C.7]) can be applied to the equation OkE^x, k) = in a small neighborhood of any 
fixed (xojKxo) £ (M.'^,'H). The unique solution is a map x ^ k^^ of class C'~^ from a 
neighborhood of Xq G M'^ to Ti. Since the corresponding geodesic Ixis) = sx + k^, 
by our uniqueness assumption coincides with 7^,. we deduce that Whence also 

X — )■ is of class C''^^. 

Re [uj) . Clearly (13.11) is a consequence of Lemma 12.11 and the uniqueness of geodesies. 
It follows from[i]) and (13. ip that S"^ is of class C'^^ on R"'. In particular, since S{x) > 
for X G M'^ \ {0}, S is of class on W'- \ {0}. If we write 7 instead of ^y^ in (13. ip . then 
using (12. 9p . (I2.10p . and the fact that d^y/dxk = se^ + dn/Oxk, with k eH, we have 

2SW^= /' (2|iG7 + 7VG.|l,)d. 

OXk Jo V (^Xk OXk J 

1 

(26^6*7 + s'jVG ■ ek'y)ds 



Thus 



1 

((2sefeG7)- + s(7VG' ■ 6^7 - (2e,.G7)-) ds 
2efcG(x)7(l). 

G(x)7(l) = S{x)VS{x), (3.4) 

and using (12. 3p we obtain 

S\x) = 7(1)G(7(1))7(1) = S\x)VS{x)G{xr'VS{x), 

from which (II. 9p follows. 

It remains to show that S is of class C' on R'^ \ {0}: From (13. 4 p we obtain 

VS{x) = S{x)-^G{x){x + k^{l)). (3.5) 

Whence it suffices to show that k^i^l) is of class C'^^. For that let us note the repre- 
sentation 

"1 , /•! 



k(1) = 2 J i^k{s) + J k{a)dajds. (3.6) 

The right hand side of (13. 6p is indeed of class C'^^. Note that in fact it follows from 
(EH) that 



d'^k{l) = 2 (d^Ks) + d^k{a)da^ds for |a| < / 



(3.7) 



where the quantities d^k and d^k are well-defined [L^{0, 1))'^— valued function, cf. H]) 
and f l212|) . 



Re Imll . It suffices to show the bounds for |x| > 1. We shall proceed by induction in 
Note that (13. 2 p for |a| = only follows from Lemma [2.11 and the representation 
k{s) = 7(s) — sx. Suppose we know the bounds for \a\ < n — 1 then we need to show 
these for |a| = n. So let a with |a| = n be given. 
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By repeated different at ion of the defining equation {d^Elx, n),h) = (for any h eT-L) 
we obtain that 

-{dlE{x, K)d''K, h) 
is a sum of terms each one either of the form 

{d%E{x, K),h) = {didl-^'E{x, K)-h)- Id = n, A; = 0, (3.8a) 
or of the form, for /c = 1, . . . , n, 

{didl+^E{x, k)- d^'K,..., d^^K, h)- (3.8b) 

k 

Id + ^ = |a| = n and 1 < <n-l. 
i=i 

Due to (11. 7p it suffices to bound the expressions in (I3.8ap and ( ]3.8bl) as 

I---I < (3.9) 

From (12. 7p we see that 

(9^9^+'=E(x,/€);/ii,...,/ifc,/ife+i) 

is an integral of a sum of {k + l)-tensors in (71 (s), ... , gk+i{s) where for each j < k + 1 
either gj{s) is a component of hj{s) or gj{s) is a component of hj{s). At most two 
factors have a "dot superscript", and if we include factors of components of 7(5) and 
factors of components of dxi-^{sx} = e,, z = 1, . . . , the total number of such factors 
is for all tensors exactly two. Whence we are motivated to group the terms into three 
types that will be considered separately below: 

A) There are no factors of components of 7('S)- 

B) There is one factor of a component of 7(5). 

C) There are two factors of components of 7(5). 

Clearly these tensors involve factors of components of d'^G{'~f{s)) also. These are 
estimated as 

l^^^^jl < C|sx|-""^'"^(l''l'^+l''l/2). ^ ^ [0, 1], (3.10) 

due to Lemma 12.11 The singular power of s in (I3.10p (depending on the a at our 
disposal) needs to be factorized into factors some of which need to be distributed to 
factors of components of hj{s) (if such components appear) and then "removed" either 
by the Hardy inequality (removing a factor s~^) 

[\-Xh{s)\^s<A\\hf, (3.11a) 
Jo 

or by the estimate (removing a factor s~^/^) 

\Hs)\ < ,/^\\h\\. (3.11b) 

Another factor of the singular power of s in f l3.10p combines with factors of components 
of dxXlis) - k{s)) = sci, i = 1, ... ,d. 

To summarize we need to look at the expressions (I3.8ap and fl3.8bp . Let hj = 
d^^ K, for j < k and /i^+i = h. After doing a complete expansion into terms (using 
the product rule for differentiation) we need to bound each resulting expression say 
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(F^^fe; hi, ... , hk+i). Recall from the above discussion that any such term is an integral 
of a (/c + l)-tensorial expression; the j'th factor is either a component of hj{s) or a 
component of hj{s). The treatment of these terms is divided into various cases. For 
simplicity we assume below that k > 1. For = we can argue similarly (although the 
treatment for = is simpler). 

Case A): Notice that we need \ri\ = \( \ + k — 1 in (13.101) . We distinguish between the 
following cases: Ai) There occur a component of hi{s) and a component of hj{s) (for 
some i j). Ail) Exactly one factor of component of hj{s) occurs. Aiii) There is no 
factor of component of hj{s). 

Case Ai): d''G{-f{s)) = s-^^^d^d^l^^^^Gisx + k{s)); \uj\ = k - 1. We choose in ( 13101) 
cr G [0, 1] such that with the given value of jr^l 

amm{\r]\, 1 + |r/|/2) = \ri\/2 =: K. (3.12) 

Upon using the pointwise bound ( I3.11bl) for k — 1 factors, the pointwise estimate 

and the Cauchy Schwarz inequality we obtain the bound 

fc+i 

\{F^,k]hi,...,hk+i)\<C\x\-''\[\\hr4. (3.13) 

m=l 

By the induction hypothesis 

m=l 

which together with (13.131) yields 

|(Fc,fc;/ii,...,/ifc+i)| < C(a;)-^(x)^/2-("-l^l)/l/i|| =C(x)V2-"/2||/,||. (g^^^^ 

Case Aii): a''G(7(s)) = s-^'^'^di'&^^^^^Gisx + K{s)y, |Ci| = |C| - 1, 1^1 = k. We choose 
a as in (I3.12p . Upon using the bound (13.1 lap for one factor and the pointwise bound 
(]3.11bl) for k — 1 factors we proceed as in the previous case using now that 

Case Aiii): d''G{-f{s)) = s-\^^\d^,^d^^^)G{sx + k(s)); IC2I = |CI - 2, \uj\ = k + 1. We 
choose a as in (I3.12p . Upon using the bound (13.1 lap for two factors and the pointwise 
bound (13.11bp for A; — 1 factors we precced as in the first case using now that 

Case B): We need |?7| = |C| + in (I3.10p . We distinguish between the following 
cases: Bi) Exactly one factor of component of hj{s) occurs. Bii) There is no factor of 
component of hj{s). 
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Case Bi): 9''G'(7(s)) = s-l'^l(9|9^(^)G(sx + k{s)); \uj\ = k. We choose in (I3J0|) 
cr G [0, 1] such that with the given value of |?7| 

amindr^l, 1 + \ri\/2) = 1/2 + \ri\/2 =: K. (3.15) 

Upon using the bound (13.1 lap for one factor and the pointwise bound (]3.11bp for — 1 
factors, the pointwise estimate 

S-^slCI,l+(fc-l)/2 < 1 

and the Cauchy Schwarz inequahty we obtain the bound 

fc+i 

h,..., hk+i)\ < ClxlVr^ J] \\hml (3.16) 

m=l 

(The factor comes from bounding |7(s)| < C\x\, cf. Lemma [2.11 ) By the induction 
hypothesis 

yielding 

|(Fc,fc; h,..., hk+,)\ < = (3.17) 

Case Bii): d^Gi^fis)) = s-l^il9^i9^(,)G(sx + /t(s)); |Ci| = |C| - 1, H = k + 1. We 
choose a as in fl3.15p . Upon using the bound (13.1 lap for two factors and the pointwise 
bound ( I3.11bp for k — 1 factors we proceed as in case Bi) using now that 

s-KsKi\s^+(^-^y^ < 1. 

For the case C) we need fl3.10p with \ri\ = \(\ + k + 1. 

Case C): a''G(7(s)) = s-^'^^di&^(^^^G{sx + /t(s)); \uj\ = k + 1. We choose in fl3:T0D 
a G [0, 1] such that with the given value of |?7| 

amin(|r/|, 1 + |r/|/2) = 1 + |r/|/2 =: K. (3.18) 

Upon using the bound (13.1 lap for two factors and the pointwise bound ( ]3.11bp for A; — 1 
factors, the pointwise estimate 

and the Cauchy Schwarz inequality we obtain the bound 

k+l 

\{F^,k;hi,...,hk+i)\<C\x\'\x\~''l[\\hml (3.19) 

m=l 

By the induction hypothesis 

k 

m=l 

yielding 
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Re[iv]). The proof is by induction in |a| and based on ( I3.5p . (13. 7p . f l2.12p andEIll. We 
notice that bound 

< for all \a\<l- 1. (3.21) 

Indeed by repeated differentiation of (12.121) it suffices, due to (13. 7p . to bound the 
L^(l/2, l)-norm of quantities of the form 

\r]\ = k + 1, > 1, ^ \f3j\ + ^ l^il = |«|- 

j=l,2 i=l,...,fc 

Using the bounds (I3.10p and [m]) we obtain that 

||5"^|Ui(l/2,l) < C(x)l/2-min(|a|-l/2,H/2)^ 

yielding f lX^ . 

Now we apply the product rule to (13.51) noticing that also derivatives of the last factor 
7(1) = x + kx{l) obey the bounds (13.211) . Using the fact that n — )• min(n — 2, (n — l)/2) 
is concave and the induction hypothesis it follows that 

Since also 

indeed the product rule (and a little bookkeeping effort) completes the induction argu- 
ment. 

Re |vj) . This is obvious from the above proofs. 

□ 

4. Proof of Theorem 11.41 111) 
We embark on proving the first assertion of Theorem 11.41 

4.1. Unperturbed case. Let Oi be the subset of order zero metrics obeying Definition 
ll.3| [T]) and let U be given as in Definition 11.11 

Lemma 4.1. Oi CU. 

Proof. Let G E Oi he given. 

Re Definition ll.ll fTi). A short calculation, using (I1.12p . gives that for all h eM.^ 

VG ■h = VP ■h + PiVG ■ hP^ + VP^ ■ hGP^ + P±GVP± ■ h; (4.1a) 

\/p.h=\^)(u\ + \u)(^\ and VP±- h = -\/P ■ h. (4.1b) 

\x\ \x\ 



Now, using (14. lap and (14.1bp it can easily be verified that 7(5) = sv is the solution to 
(12.121) that satisfies 7(0) = 0, and 7(0) = v, for any given v G M*^. By standard unique- 
ness of the solution to an initial value problem for an ordinary differential equation this 
7 is the unique solution to (12.120 that satisfies 7(0) = 0, and 7(0) = v. Since E £x is 
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a geodesic of G if and only if 7 satisfies (]2.12p . 7(0) = 0, and 7(1) = x, then 7(3) = sx 
for s G [0, 1], and thus G satisfies Definition II .ll fTl). 

Re Definition ll.l| [2ll. We will use the representation ( ll.Sp . Since for all x G M*^ the 

geodesic is 7(5) = sx, a not very short but elementary calculation, using (14. lap . ( I4.1bl) . 
and the shorthand notation 7 = d'j/ds, h± = P±h and 7 ■ a; = (7, cj), yields 

4 / ^VG-hhds 

4/ {h^ ■ - h^Gh^}p\ds (4.2) 







and 



1 

2 



7(V'G; h, h)'yds 



Jo I I7I I7I J 



1 U,|2 



= 2/ {-\h^\^ + h±GhA[^ds. (4.3) 
Jo I7I 

Thus, using for the first and last equations below that 

■ h^m^i - 2|/.j2|7r/i7r = mimm, (4.4) 

we find 

{dlEix,0)h,h) 

= 2 [ [hGh - 2hxGhX^^ + h^GhM-Xds (integrating by parts) 
Jo I I7I 

= 2 [ {\Ph\' + s'{s-'h^yG{s-'h^y}ds (using I7I/I7I = s"') 
Jo 

>2 [ {\Ph\^ + as^\{s-^hj_y\^}ds (hj 
Jo 

= 2 (|P/ip + a\{h±y\^)ds (integrating by parts). (4.5) 
Jo 

We conclude (using the bound a < 1) that 

{dlE{x,0)h,h) >c\\hf; c<2aandhe'H. (4.6) 
Hence G satisfies Definition II. 1| [2|) and therefore Oi C W. □ 
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We will extend the above proof to the case of geodesies 7 G for metrics "near O" 
a priori not knowing that geodesies are unique. To do this we need dynamical control 
of the geodesies and this will be provided under the additional condition fll.lSp . First 
we discuss the case of G E Oi (as in Lemma 14.11) . Let 7 denote any non-constant 
(maximal) geodesic for such a metric (note that we are here not assuming 7(0) = but 
only the differential equation fl2.12p ). Introduce the observables 

A = ■ 7 and 5 = 7-7; (4-7) 

niG 



Note that 



A = q-T. g= ^'|'[^^^ , (4.8) 
|7|g 

cf. (11.121) . In particular 

J^,B^<\. (4.9) 

Lemma 4.2. Suppose G E Oi and that 'y is a corresponding non-constant geodesic. 
Then 



A=^^(Px(7)7,TPx(7)7); T = G(7) + 2-S-VG(7). (4.10) 

I7I |7|g 

In particular if G E O 



A>c'-j^[l- A'y (4.11) 
with c > given by (I1.13P . 

Proof. Using fl2.12p . ( 14. lap and (]4.1bp we compute fl4.10p . Note that the denominator 
171(3 is preserved, cf. fl2.13p . 
As for (I4TT]) we use fgrTOj) . f lLTS]) and 

\pm? = 1 - \p±m\' = B' = (4.12) 

Lemma 4.3. Let G G Oi, and let A and B he given by fji??]) for any (7, 7) G (]R"'\{0})^. 
We have 

b{l - B^) > 1 - (4.13a) 

a-\l- A^)>1- B^. (4.13b) 

Proof. Using f ll.l2p and f l4.12p we can estimate 

^ = iPi'p + \P^G'/^P^^\' < 1 + (6 - 1)(1 - P2), (4.14) 

nr 

which in turn using (14. 9p . fl4.12p and the fact that b > 1 yields 

l-B^>l-A'^-{b-l){l-B'^). (4.15) 
Obviously (I4.13ap follows from fl4.15p . 
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Next we mimic the proof of ( I4.13ap . We have 1 — = 1 — ||p--B^, and letting 
q = G^^'^'j/\'j\G (as in fl4.8p ) we estimate 

I ■ |2 

n\ = |Pg|2 + |P^G-V2^|2 < 1 + (a-1 - i)(i - A^), 

niG 

cf. (14.141) . Whence the analogue of (I4.15P holds, and we conclude ( I4.13b[l . □ 

We remark that only fl4.13bl) will be needed. The estimate fl4.13ap is given only for 
completeness of presentation. 

4.2. Perturbed case. Now let G G O be given. We shall use the notation to 
denote any metric of order zero obeing — < e. The positive parameter e is an 
order parameter which we will take sufficiently small, say e < eo, in terms of quantities 
given by the fixed unperturbed G. We shall use the observables A and B of (14. 7p 
defined in terms of G but now evaluated at 7 — > 7^; clearly they are well-defined for 
7.(s),7e(s)^0. 

Lemma 4.4. There exist eo > and Ci,C2,C^ > such that if \\G^ — G\\i < e < eo 
and 7e is any non-constant geodesic for the metric G^ with 7^(0) = 0, then 

i) 7,(s) 7^ for all s > 0. 

ii) B = 5(7,(s),7,(s)) >l-eCi for all s > 0. 
Hi) There are bounds 

(1 - eC,)\%m' < Ms)\' < (1 + eC2)|7.(0)r, (4.16a) 

(1 - eC,)\s%m' < |7.(s)r < (1 + eC,)\s%m\ (4.16b) 

Proof. Re [1]). Let 7^ be any such geodesic. Being non-constant implies 7e(0) 7^ 0. 
Indeed note at this point that the quantity 7eGe(7e)7e is constant. Moreover, using also 
fO]) and the fact that G(0) = /, 

a|7e(3)P <|7.(5) <&|7.(^)P, (4.17a) 

(1 - eC)|7.(0)P < |7.(^)l^(,.(.)) < (1 + eC)|7e(0)r, (4.17b) 

(1 - eC)b-'%{0)\' < Ms)\' < (1 + eC)a-'\%m'. (4.17c) 

We shall use fl4.17ap and fl4.17bp later in the proof while (14.1 7cp is stated just for 
completeness (note that fl4.16ap is stronger than fl4.17cp ). 

Whence the observables A and B are well-defined at this geodesic on an interval 
of the form ]0,so[. We need to show that so can be taken arbitrarily large. Suppose 
not, then we let sq be the first positive nullpoint, 7e(so) = 0, and we need to find a 
contradiction. For that it suffices show the bound 

A{s) >l-eC for all s g]0, So[, (4.18) 

where the constant C > depends only on G (i.e. it is independent of e, G^ and 7^). 
In particular we can assume that eC < 1. It then follows that A,B>0. Whence by 
the computation 

f|7,|2 = 2(7„7,) = 2|7,||7,|i?>0, (4.19) 
we see that s — )■ |7e(s)| is increasing yielding the contradiction and whence showing [1]). 
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It remains to show (I4.18p . We compute the time- derivative of A and find the following 
extension of (14. lip . 

A > c^-^^\^^^ (1 - A') - ec\^. (4.20) 

l7e| 17^1 

Note that 



and whence due to fl4.17ap that 



-1 



The constant C contributes to the constant C in fl4.20p . Another term of the form 
of the last term in f l4.20p comes from comparing the right hand side of the geodesic 
equation fl2.12p for % with the same expression replacing G'e(7e) — )■ G(7e). 
Using (I4.17ap again we can simplify ( I4.20p as 



A > -cV^nA {l-A^-eKy,K=^^. (4.21) 



Combining f l4.2ip and the fact that lims^o+ ^{s) = 1 (here once more using that 
G(0) = /), we obtain f l4.18p for any C > K/2 provided e is chosen small enough 
(precisely we need 2C — K > eC"^). 

Re [u]). We write 5 = 1 - (1 - 5) and then use fl4.13bl) and f liTTSD (we use fHTISD with 

So = oo) to estimate 

(1 + B)B > (1 + 5) - tCa-^{l + A)>l + B- aCa-^. (4.22) 
We subtract B from f l4.22p and obtain 



> V 1 - e2Ca-i > 1 - tCa-^ -e^C>\- eCi; Ci := 2Ca~\ 

Re[iii]). Using [n]) we obtain, cf. f Hlij) and fH:22|) . 

(l-62C0l7e(^)r<|7.(^)l^(,.(.))- 
Also, using ( I4.18p . we have 

(1 - eC)|7e(s)|G(7.(s)) < l7e(s)|G(7.(s))^ < l7.(s)l, 

yielding 

l7e(^)|^(,.(.))<(l-eC')-2|7e(.)p. 

Using the shown two-sided estimates in combination with (]4.17bp we obtain fl4.16ap . 
Using 

f l7.l = 

cf. fl4.19p . we obtain ( I4.16bp from f l4.16ap andEI]) by integration. □ 
We have the following version of Definition I1.1| |2|) for the perturbed metrics. 
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Lemma 4.5. There exists eo > (possibly smaller than the eo of Lemma 4_;4_) such that 
if \\Ge — G\\i < e < eo, x G M'^ and 7^ is any geodesic for the metric G^ emanating from 
with value x at time one, then writing 7e(s) = sx + k{s), denoting by the energy 
functional (11 .4^ with G replaced by G^ and using the positive number a given by ( \1.3\\ 
for the metric G, 

{dlE,{x,K,)h,h) >a\\h\\^, heU. (4.23) 

Proof. Note that for a; = the geodesic 7^ is unique due to Lemma 1131 E]), it is given 
by 7, = 0. By ([L8l) ^^^^(0, 0) = 2G'e(0), so obviously Km holds in this case provided 
a < 2a-2ed (< 2G,{0)). 

Suppose next that x ^ 0. We can mimic the proof of Lemma HIT] using Lemma 1^741 By 
letting s = 1 in (14.16b|) we obtain from (I4.16ap and ( ]4.16b|) that there exist C4, C5 > 
such that 

(1 - eC4)|xp < |7,(s)|2 < (1 + eC4)|x|2, (4.24a) 

(1 - eG,)\sx\^ < |7.(s)|' < (1 + eG,)\sx\^. (4.24b) 

A consequence of fl4.24ap and (]4.24bp is that the quantity |7e|/|7e| appearing when we 
try to repeat the proof of Lemma [4.11 effectively is given by to be used in the last 
part of the proof only. More precisely 

(1 - eGe)s-' < |7,|/|7.| < (1 + eCe)s~\ (4.25) 

We calculate using (14.1bp . fl4.12p . Lemma I^^ HIjl and the notation u = 7e/|7e| 

P±7. = 0(v^)|7.|, (4.26a) 



ds 



Pile 



{UJ\ + \UJ) 



%\ / \ 17, 



0{V~e) g^; (4.26b) 



here and below the notation 0{^/e) is used for any (matrix- valued) function of s G [0, 1] 
obeying \0{y/e)ij\ < y/eG uniformly in s, x, G^ and 7^. In fact we can above choose 
G = 72 ^ and G = 2^f2n[, respectively. 
From (I4.26bp we obtain 

I (P(7e)/^) = P{n.)h + O(v^) M/,. (4.27) 

Below we stick to the notation = (/i_l) , although by fl4.27p we could have choosen 
an alternative interpretation. 

Now, using fHlap . (l4Tb|) . ( I4.26ap . fOTj) and Lemma |43] In]) the analogue of fOj) 
reads 

4 I 7,VG, ■ hhds 



= 4/ {/ix -/ix -/i±G'/ix + /iO(v^)/i + /iO(x/^)/ip^}P^ds. (4.28) 
Similarly the analogue of (14. 3p reads 

\ %(V^G,;h,h)%ds = 2 [ { -\h^\^ + h^Gh^ + hO(^/^)h]\^ds. (4.29) 

Jo Jo Hel 
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Next we compute using (I2.12p . (14. lap . ( 14.1bl) and Lemma lOlEI]) 

(l7.|/l7.|)- = -(l7.|/l7.|)^(l + 0(v/i)). 
Whence the analogue of (14. 4p reads 

■ /^±|7e|/|7.| - 2|/.j2(|7e|/|7e|)' = 2 (/^l |7. | / | )■ + /^O ( v^i) ^ 1% | / |7e | (4.30) 

We insert (I4.28P and (I4.29P into (II. 7p . integrate by parts using (I4.30p and obtain the 
following partial analogue of (14. 5 p 

{dlE,{x,H)h,h) (4.31) 

= 2 I |/iG'e/i-2/ixG/ixN + ^±G'^±W + ^C»(v^)/i^ + /iO(v^)/i^lc/s. 

Finally we invoke (I4.25P and (I3.11ap . and obtain from (I4.3ip using again (I4.27P the 
following analogue of (14. 6 p 

{dlE,{x,K)h,h)>2 [ {a- V^C)\h\^ >a\\hf. (4.32) 
Jo 

□ 

Proof of Theore'm \1.4\ \^- Let $e = expQ^(l-) : TMq — )■ M'^ denote the exponential map 
for the perturbed metric (close to a given G G C) at the point G M"^ and evaluated 
at time one. The positivity of the Hessian along any perturbed geodesic as guaranteed 
by Lemma [4.51 implies that $e is a local diffeomorphism, cf. |Chl Theorem 2.16] and 
pn Theorem 14.1]. By Lemma O $e maps TM^ onto W^. By (I4.16bp is proper, 
whence $e is one-to-one, cf. |Be| Theorem 5.1.4]. We have verified the uniqueness 
property of Definition ll.ll fTj). 

The uniform positivity property of Definition ll.l| [2|) follows from Lemma 14.51 (with 
c = a in (II. 7p ). Whence indeed G W for e < eo. □ 

5. Proof of Theorem ll.4l ini) 

We shall prove the bounds (I1.15ap . (I1.16P and (I1.15cp . Let G G O be given. Using 
the convention of Subsection 14.21 we write G = for perturbed metrics, and we shall 
again require \\G^ — G\\i < e < ior some small eo > 0. 

Re (ll.lSap . By the variational definition (O]) and (ES]), ^^(x) - |xp = 0(e)|xp 
uniformly in a; G M*^ \ {0}. Whence we obtain uniformly in x G M'^ \ {0} 

S{x) - \x\ = {S\x) - \x\^)/{S{x) + |x|) = 0(e)|x|, (5.1) 

yielding (I1.15ap by division with 

Re (fTTeb . We compute 

d^s{x) = - (5.2) 

\x\ \x\'^ 

For the first term we use (13.51) , (I1.15al) , (I4.24al) and Lemma I4.4| [iip to write 

\x\d,S{x) =Xi + 0{,/^)\x\. (5.3) 
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As for second term we use (ll.lSap to write 

^^ = a;. + 0(e)|a;|. (5.4) 

Clearly ffTTel) follows from (Q-dMl). 
Re ( ll.lScj) . 

Step I. We prove the uniform bound 

\\d"K\\ < v^C for |a| = 1. (5.5) 

Note that (15. 5p without the factor ^/e to the right follows from Proposition 13.11 Due 
to Lemma 14.51 it suffices to bound the expression (13. Sal) for n = 1 as 

\{d^d^E,{x,K),h)\ < yAC\\h\\ for a = ej. (5.6) 

As in the proof of Proposition 13. II we compute the x-derivative using (12. 7p yielding four 
terms. Up to an error of order 0(e)||/i|| we can replace — )■ G, and whence it suffices 
to show that the sum of the following four expressions is of order 0(A/e)||/;.||: 

Ti = / 2ejVG ■ h% ds, 



T2 = 2hGej ds, 
Jo 

T3 = s'jfVdjG ■ h'j^ds, 
Jo 

Ti = 2shdjG% ds. 
To do this we use fl3.11ap . I KTa^ . (liTbl) . fl4A2|) . ( OSj) and Lemma gain]) and obtain 

T^= [ 2e,{P^-P^GP^)^Js + 0{V~^)\\h\\, 
Jo 

T2= [ 2ej{P + P±GP±)hds, 
Jo 

n = [ 2e,{P^GP^~P^)i^ds + 0{V^)\\hl 
Jo 

T4= /" 2ej{P±- P±GP±)hds + 0{y/^)\\h\\. 
Jo 

Clearly it follows that 

Ti+T3 = 0(v^)||/i||, 

T2 + T4= / 2ejhds + 0iy^e)\\h\\=0iV~e)\\h\\. 
Jo 

Whence we have proved (15.60 . 
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Step II. We shall prove the uniform bound 

\d:f^il)\ = \d:U^) - e,| < VeC for a = e,. (5.7) 

We claim that 

djk{s) = s-^F{s) where [ \F{s)\^ds < eC\ (5.8) 

Jo 

Note that (15. 7p follows from (13. 7p . (15. 5p and (15. 8p . The difficulty here is not to show 

that the quantity sdjk{s) is square integrable but rather to show that its L^-norm is 
bounded by y/eC as stated in (15.80 . We shall show that 

djk = G-^PiVG ■ %P± {s-^djK - djk) + s-^F where \\F\\l2 < v^C, (5.9) 

which combined with (I4.25P and (15. 5 p supplies (15. 8p . 

To prove (15. 9 p we proceed as follows. Differentiating (I2.12p with respect to Xj and 
using (I4.25p . (I3.11ap . (15. 5 p and the facts that (9^7^ = scj + djti and \\G^ — G\\i < e we 
obtain 

G.djk = dj%VG ■ e.7, + 2-^,djVG ■ e,% - djVG ■ 7,7e (5.10) 
- VG ■ dj%% - VG ■ %d,% - djGk + 0(e)/s, 

where here and henceforth 0{eP) stands for a function with L^-norm bounded by e^G. 

In the remaining of the proof of ( 15. 9p we will repeatedly use (14. lap . (14.1bp . (14. 250 . 
(gSSil), flXTTil) and ([53]). 

First we estimate djGk. To this end we note as above that 

djGk = djGG-\2-^%VG ■ e,% - VG ■ %%) + 0{e)/s. 

We compute 

djGG~^,VG ■ e.7, = VG ■ dj%G-^[y,VP ■ e.7, + 7eP±VG ■ e,P±% 

+ 7eVPx ■ e.GP^7. + ieP±GVP± ■ e.7.] 
= VG ■ dj%G-'%VP ■ e.% + 0{V~e)/s 
= 0{V~e)/s. 

Similarly we obtain djGG^^VG ■ 7^7^ = 0{^/e)/s, and therefore 

djGk = 0{V^)/s. (5.11) 
Now we estimate the first five terms on the right side of (]5.10p . 
(i) For the first term we have 

dj%VG ■ e.% = dj%[VP ■ e. + P^VG ■ e.P^ + VP± ■ e.GP^ + P^G\/P^ ■ e.]7, 
= dj%WP ■ e.7. + d,%P^GWP^ ■ e.7e + 0{^fe)/s 
= ((7„7,)/|7,|)((9,-7e,P±e.) - {GP^d,%,P^e,))+0{^e)/s. 

Thus, 

a,-7.VG ■ e.7e = ((7e,7e)/|7e|)(/ - PLG)P^d,% + 0(v^)/s. (5.12) 
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(ii) To estimate the second term we consider 

%djVG ■ e,% = 7.V[VP • e. + P^VG ■ e,P^ 

+ VP± • e,GPi_ + P±GVPi_ ■ e.] ■ 5,7,7, 
= 7.V(VP ■ e.) ■ 9,7.7e + 7eV(VPxV ■ e.GP^) ■ d,^,% 
+ 7.V(PxGVPx ■ e.) ■ 9,7.7. + 0(v^)/s. 

It follows that 

%djVG ■ e,% = 7eV(VP ■ e.) ■ 8^-/,% + 7, VP ■ e.GVP ■ 9,-7e7. 

+ 7.VP ■ dj%GVP ■ e.% + 0{^t)/s. (5.13) 

Using the fact that for any fixed vectors h and z we have V(x) ■ h = P±h/\'j^\ and 



P±z 



\x\ 



1 



[{x, z}P±h + {P±h, z)x + (x, h}P±z] , 



we find that 

7,V(VP ■ e.) ■ 5,-7,7, = 7,V 



(7e| + I7.) 



2a, 7.) (7., V (^^) ■ 9,7e) + 0{yre)/s 



-2(7., 7.) (7 



{P±dj^e, e.)7. 



+ 0(v^)/s. 



Hence 



7,V(VP • e.) ■ 9,7e7. = -2%^P^9,7. + OiV~e)/s. 



Moreover, using twice fl4.1bp we obtain 

P±5,7, 



7.VP ■ d,%GVP ■ e,% = % 



{%\ + lie) 



GVP ■ e.7. 



a, 7.) ( ^^^4^, VP ■ e.7. > + 0{V-e)/s 



= {{%,%)V\%\')P±GP^d,j, + 0{V~e)/s. 

A similar calculation supplies 

7. VP ■ e.GVP ■ 9,7.7e = ((7e,7e)Vl7.r)^±G'Px«9,7e + 0{V~e)/s, 
and putting the last three estimates together in fl5.13l) gives 

2-^,d,VG ■ e.% = ((7„7.)Vl7er)(P±G - I)P±d,^, + 0{V~e)/s (5.14) 
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(iii) Concerning the third term of (I5.10p we have 

d,VG ■ 7e7e = V(VP ■ 7e) " 5,7.7. + V(PxVG ■ %P^) ■ d,%% 

+ V(VPx ■ %GP^) ■ 57.7. + V(PxG'VPx ■ 7.) • d,%% 
= V(VP ■ 7e) ■ 9,7.7. + ^±VG ■ 7.VPx ■ 9,-7.7. 
+ P^GViVP^ ■ 7,) ■ 9,7.7. + 0(v^)/s. 
Analogous calculations to the ones leading to (15.141) yield 

V(VP ■ 7.) ■ 9,7.7. = -((7.,7.)/l7.|)'^±5,7. + 0(v^)/s, 
P^VG ■ 7.VPx ■ 9,7,7, = -((7.,7.)/l7.|)^±VG • 7.Px9,7, + 0(v^)/s, 

and 

PxGV(VPx ■ 7.) ■ 9,7,7, = ((7.,7.)/|7.|)'^±G'Px5.-7. + 0{V^)/s. 
Therefore 

9,-VG ■ 7.7. = - ((7.,7.)/l7.|)'^±fi',7. - ((7., 7.)/l7.|)i^±VG • 7.P±9,7. 

+ {{%,%) /\^,\yP^GP^d,% + 0{V~e)/s, (5.15) 

(iv) To estimate the fourth term of (I5.10p we proceed as follows 

VG ■ 9,7,7. = (VP • 9,7. + P± VG ■ 9,7,Px 

+ VPx ■ 9,7,GPx + PxGVPx ■ 9,7.)7. 
= VP • 9,-7, + P^GVPx ■ 9,7.7. + 0{V~e)/s 
= ((7.,7.)/l7.|)(/-i^±G)Px9,7. + 0(v^)/.. (5.16) 

(v) For the fifth term of (15.101) we have 

VG ■ 7,9,7, = (VP ■ 7, + P±VG ■ 7,Px 

+ VPx ■ %GP^ + P^GVP^ ■ 7,)9,7. 
= PxVG ■ 7.Px9,7, + 0{^t)/s. (5.17) 
Substituting fISlTD . fl5:T2|) . (jEHD, fl5J[5|) . f lSTTG]) and f lSTTj) in f lSTTOj) we find that 

G,9,^ = P^VG ■ 7,Px (((7„7,)/|7.|)5,7. - '9,7,) + 0{^e)/s, 
and using now ( 14.25^ and Lemma [4.4| [iil) we arrive at (15. 9p due to cancellation. 

Step III. We shall complete the proof of (I1.15cp by differentiating (15. 2p to obtain a 
representation of the second order derivatives and then using (I1.15ap . (I1.16p . (13. 5 p and 
fl5TD . We have 

Pi — ^ii^ -^j^i^ + XjdjS 6ijS ^ ^ XjXjS ('cloA 
OijS[x) - —— — TTF TT^' ^ ^ ^ 

Using f lTTSal) and f lLTej) it follows from fl518|) that 



9i,S' XiX,- - + 0(v/e) 
9i,s(x) - -p^ = — ^ f— . (5.19) 

\x\ \xr 
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As for the first term on the right hand side of (15.180 we use ( I3.5P and (15.70 to write 

_ S-'G{e,+0{^e)) ^ {d,S-')G^,{l) ^ g'^VG ■ e,7.(l) _ 

The ith component of the right hand side of (I5.20p is of the form 
Obviously 

I 7' I I nf* I ^ I I ^ 

By adding these expressions we obtain that 

dijS ^ XiXj — 6ij + ©(-y/e) 2 

We combine (I5.19P and ( I5.2ip and conclude (I1.15cp . □ 

Remark 5.1. An analogue of (15. 5p was proved by a different technique in jBa] (see |Bat 
Lemma 3]) which does not work in the present context. Note that in [Ba] also an 
equation like ( 15. 9 p is used, and in fact it is used to prove the analogue of (15. 5p . Howewer 
in our case there is no smallness of the factor G~^PxVG ■ 'jeP± and consequently the 
technique in [Baj is not apphcable. Note also that |Bat Lemma 4] is an analogue of 



6. Proof of Theorem 11.61 

In this section we need the Sobolev spaces T-L^ := PFo^'^(0, 1)"^, 1 < p < oo, consisting 
of absolutely continuous functions : [0, 1] — )■ M'' vanishing at the endpoints and having 
h G -^^^(0, 1)*^ = L^QO, 1[, M*^) (throughout this section we use the notation for this 
vector- valued space). The space "H^ is equipped with the norm 

\\h\\nr = \\h\\p= (^j\h{s)\Pdsy''. (6.1) 

The first goal is to find a substitute for the positivity bound (I4.23p used in the proof 
of Theorem 11.41 We shall use the observation, of. the proof of Lemma 14. H that 

{dlE{x,0)h,h2) = [ 2s\{s-^hi{s)yG{sx){s-^h2{s)y}ds. (6.2) 
Jo 



Motivated by (16. 2p we develop some functional analysis which then will be applied to 
an extension/modification for perturbed geodesies. 
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6.1. Hardy type bounds and duality theory. 
Lemma 6.1. For all p g]1, oo[ there exists Ap > such that 

\\h{-)/s\\p =[j^ \h{s)/s\^dsy'' < Ap\\h\\nP for all h G 71^. (6.3) 

Proof. We refer to |HLPl Theorem 327]; the bound is valid for Ap = p/{p — 1). □ 
Lemma 6.2. For all p g]1, oo[ there exists Bp > such that 

\\h\\w < Bp\\hi-) - hi-)/s\\p for all h G U^. (6.4) 
Proof. Consider the linear maps 

LP{0, 1^3 f^ Sf; {Sf){t) = t-' f f{s)ds. (6.5a) 

LP{0,lY3f^Tf; {Tf){t) = -t s-^f{s)ds. (6.5b) 

Note that 5* is bounded on with ||5'||e(LP) < Ap, cf. Lemma l6.ll and its proof. 
Next note 

|(r/)(t) = ri(r/)(t) + /(t), (6.6) 

and recall the standard fact 

L^(/, g) = (L^(/, gy if q-' = l- p-'- (6.7) 

in our case / =]0, 1[ and g is the Hilbert space W^. In terms of ( 16. 7p we can rewrite 
dnH) as 

f:=jfr = -S* + 1, 

and since S G B{L'^) we conclude that T G B{Lp). 

Finally noting that for any h G "H^ we can write h = Tf where f{s) = h{s) — h{s)/s 
we obtain (lO) with Bp = Ag + 1. □ 

Remark 6.3. By the estimate (due to the Holder inequality) 

\ l\s-'fis)ds\<t-'/^ip-lY-'/^f\\p, 

it follows from the proof of Lemma (6.21 that the operator T of ( I6.5b[l is in B{U' ,1-1^) 
with norm ||T||g(iP) < Bp. In fact (using here also Lemma [6.11) T : W ^ 'W is a linear 
homeomorphism. 

Lemma 6.4. Suppose G is a given continuous function, [0, 1] — )• G{s) G iSd(M), using 
here notation of Subsection \l.I\ Suppose there are constants a,b> such that 

a\y\^ < yG{s)y < b\y\\ y e R'^ and s e [0,1]. (6.8) 

Introduce for h G "H^ and g G 'H'^, + p~^ = 1, the pairing 

[h,g] = [\s'{{s-'his)yGis){s-'g{s)y}ds, (6.9a) 

^0 
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and the associated quantity 



np,G= sup \[h,g]\. 

ll9llTi9<l 



(6.9b) 



This quantity \\ ■ ||-^p,g is a norm on and 

Ci \\h\\ (6.10) 
C, = 26(1 + Ag){l + Ap), C2 = {2a)-'BgB,. 

Here Ci and C2 are given in terms of the constants of Lemmas \6.1W6.2\ and i \6.8\\ . 

Fmally {W)* = W m the sense given by the pairing fl6.9al) . This means that for all 
I G (T-L^Y there exists a unique g G 'H'^ such that 

1(h) = [h, g] for all h G H^, (6.11) 

and vica versa any g G "H"^ defines an element I G {H^Y (16.111) . Moreover the 
identification is a linear homeomorphism. 

Proof. Notice that indeed the pairing (I6.9ap is well-defined due to the Holder and 
Minkowski inequahties and Lemma 16.11 Using these bounds we also get the first esti- 
mate of (EHOj). 

To prove the second estimate of fl6.1UI) we first use Lemma 16.21 and (16.71) to estimate 

>i 

SIS 'ft(s) I -iU[s]T[s]as 



HP < B 
< B 



p sup 

li2G/||i<,<l 



p sup 

ll/llL9<(2a)-i 



(s-i/i(s))2G(s)/(s)rf, 
s(s-^/i(s))2G(s)/(s)ds 



Next we introduce for any f & L'^ the function g = Tf which according to Remark 16.31 
is an element of "H^. We obtain 

"1 



/i o/p ^ B 



p sup 

\\f\\Li<{2a)-\g=Tf 



< Bp sup 

\\g\\^g<(2a)-^Bg 



Bp{2a) ^Bq\\h\\-Hv^G- 



s{s-^h{s))'2G{s)s{s-^g{s))'ds 
s{s-^h{s))'2G{s)s{s-^g{s)y ds 



We have proved (16.101) . 

The identification asserted next follows similarly from (16. 7p and Remark 16.31 The 
bi-continuity is a consequence of (I6.10p . □ 

Corollary 6.5. Let G be given as in Lemma 6^_, and define P G B{'H), H = 'H? , by 

{Phi,h2) = [hi,h2]; /li,/l2 GH. 

Let Ci and G2 be the constants in (I6.10p . 

For all p G [2, cxo[ the operator P G B(7i^) and it obeys 



P\\b{hp) ^ 2Ci. 



(6.12a) 
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In particular (with q g]1,2] being the conjugate exponent) 
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\{Ph,g)\ < 2Ci\\h\\nr> Mw for all heW 


and g E IH? . 


(6.12b) 


Conversely, if 






h G T-i^ and \ {Ph,g) \ < c (7 ^9 for all 




(6.13a) 


then 






\\h\\w<cC2. 




(6.13b) 


Proof. We show that Ph G W for all h eW li p e [2,oo[. 


Fix any such p 


and such h, 



and let q g]1, 2] be the conjugate exponent. Then Ph G W if the expression 

||P/i||^p = sup I / {Ph)'{s)f{s)ds\ < 00. 
We introduce the map 

L^3f^g = Rf; (i?/)(t) = 



fis) 



f{l)dl)ds. 



(6.14) 



Clearly R maps into 'H'^ and in fact (due to the Holder inequality) 



\R\ 



i3(Li,ni) 



< 2. 



Whence 



ll^^ll 



HP 



sup 

ll/lli,Q<i,/eL2 



iPhyis)iRfns)ds 



iPhyis)g{s)ds 
\[h,9]\ 



< sup 

!lsllw9<2,(;e-H2 

< sup 

Il9lk9<2,(;ew2 

<2Ci\\h\\w. 

So indeed Ph G Ti^ and the estimates show fl6.12ap . 

Clearly f l6.12b[) follows from the Holder inequality and f l6.12ap . 
As for fl6.13bp we note that by a density argument and f l6.13ap 

\\h 



\np,G= sup \[h,g]\ 

\\9\\ni<'i-,9&H2 



sup \{Ph,g)\<c. 



Whence we can conclude by f l6.10p . 



□ 



Remark. In our application, see Subsection 16. 2[ G{s) will be a composition of a metric 
in O and a perturbed geodesic (as in (16. 2p for the unperturbed case). We shall conclude 
(]6.13bp for functions in question by verifying (I6.13ap . 

For completeness of presentation let us note the following additional property of 
the operator P (shown to be in B(7i^)): For all p G [2, 00 [ in fact P is a linear 
homeomorphism on Ti^. Outline of a proof: We verify the condition fl6.13ap using the 
Holder inequality and deduce from the conclusion fl6.13bp that 

-HP < C2\\Ph\\nP for all h G H^. 
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Whence P : — Ti^ is injective and has dense range. It remains to show that its 
range is dense. For that we note that 



{h, g) = [ kg ds, heTi^ andg e Ti", 
Jo 



defines a pairing giving another identification (T-L^Y = W (use the first step of the 
proof of Corollary 16.51 and Hahn-Banach theory). Whence if the range is not dense we 
can pick g E such that {Ph, g) = for all h G "H^, violating that {Ph, g) = [h, g] 
and the last part of Lemma 16.41 

6.2. Bounds of in "H^, p G [2, oo[. We shall improve on Proposition 13 . 1 in the 
case of a metric G E M. close to a given G G C (as in Theorem II. 4liip ). We shall use the 
convention used in Subsection 14.21 and write in terms of an "order parameter" G = 
if \\G^ — G\\i < e. Following the proof of Proposition I3.1|[iii|l we can for 1 < |a| < / — 1 
represent the quantity 

-{dlE,{x,K)d'^K, h) 

as a sum of terms involving derivatives d^K where < |a| — 1. Using a new induc- 
tion hypothesis we shall estimate these terms individually improving the corresponding 
bounds of the proof of Proposition 13. lllml) . On the other hand we can write 

dlE,{x,K) =P + R; (6.15a) 
{Ph,h2) = [h,h2] := [\s\{s-'h,{s)yG{-f,{s)){s-'h2{s)y}ds. (6.15b) 







The second term in f l6.15ap can be estimated by the following 

Lemma 6.6. Let p G [2, oo[ and q be the conjugate exponent. There exists eo > and 
C > such that for all x eW'' and for ||Ge - < e < Cq 

\{Rhi, /i2)| < -\/eC||/;,i||^p 11^211^9 for all hi E Tip and /i2 G 7^2- (6.16) 

Proof. The result follows from f l4.3ip using f l4.25p . fl4.27p . the Holder inequality and 
Lemma 16.11 □ 

Lemma 6.7. Let p, q, eQ> and C > be given as in lemma l67Si Suppose 1 < |a| < 
/ — 1 and that for some constant > independent of x eM.'^ and e G [0, cq] 

\{dlE,{x,K)d''K,h) \ <C^{xy-^''^\\h\\n', for allhEn^ (6.17) 

Let a,b > be the constants from (11. 3p determined by the metric G and let C2 be the 
constant from (I6.10p (fixed in terms of a andp). Then 

(1 - v^CCa) ||9"/€||^p < C«C2(a;)^-l"l. (6.18) 

Proof. Clearly (16. 8 p holds for the example G'(7e(-)) used in fl6.15bp . From (I6.15ap . 
fl6.15bp and Lemma [6.61 we deduce (I6.13ap with 

c = C^{xY'\''\ + ^eCW^Wn., 
and we obtain from f l6.13bp that 
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yielding f l6.18p by subtraction. □ 

Proposition 6.8. Let p G [2, oo[. There exist eo > and Cp > such that for x G M*^, 
\\G, - G\\i < e < and I <\a\ < I - I 

\WnU.<Cp{xY-\"\. (6.19) 

Proof. Using fl6.18p for all sufficiently small e we only need to demonstrate f l6.17p . So 
let us verify ( I6.17P for the case | a | = 1 : We mimic Step I in the proof of Theorem ll.4liip 
using the Holder inequality and Lemma 16.11 (to replace the Cauchy Schwarz and Hardy 
inequalities, respectively). This yields f l6.17p for |a| = 1, in fact with any extra factor 
y/e. We have shown the proposition for the case 1 = 2. 

Next let us suppose / > 3 and that we know the bounds (]6.19p for all p G [2, oo[ and for 
I a I G [1, n — 1] where l — l > n >2 (notice that eo > may depend on p). Let p G [2, oo[ 
and a with |a| = be given. The proof is complete if we can show the existence of 
eo > such that (16.171) holds for the given p and a provided \\G^ — G\\i < e < Eq. For 
that we need to modify the proof of Proposition l3.1liiT|l . Again we can assume that 
\x\ > 1. We shall use the induction hypothesis for 

p ^ p := Anp. (6.20) 

This particular choice fixes some eo > that we claim indeed works for (I6.17P (it is not 
claimed to be an "optimal" choice of p although growth in n is indispensable). The 
Hardy inequality (13.1 lap is replaced by Lemma 16.11 and (I3.11bp by 

\Ks)\<s^-^'nh\\ni'-. (6.21) 
here we shall use the p of (I6.20p . We shall also need the generalized Holder inequality 

/ \fi{s)\ ■ ■ ■ \f^{s)\ ds < ll/illp, ■ ■ ■ for IM + 1/p^ < 1. (6.22) 
Jo 

Finally we shall use the following special case of (13.101) 

\d"gij \ < Clsx]-!^"!; (6.23) 

here and henceforth we omit the subscript e (obviously G is independent of e). Note 
that (I6.23P is the best possible bound at infinity. The main issue compared to the proof 
of Proposition I3.1|[iiil) is that the improved pointwise bound (I6.2ip . used to h = d'^^ k, 
compensates for the worse singularity at s = when using (I6.23p . 

Now let us look at some details: There are cases A), B) and C) defined as in the 
proof of Proposition I3.1|[iiil) (and we treat only k > 1). 

For the case A) we need (16.230 with |?7| = \(\ + k — 1, and we consider subcases 
Ai)-Aiii) defined as before. 

Case Ai): (9''G(7(s)) = s-l«l(9^(9;^(^)G(sa; + k{s)); \co\ = k - 1. Suppose ffist that 
k > 2 and i,j < k. Then we shall use (16.220 with m = 4, pi < p/n, p2 = p^ = p and 
P4 = q, fi{s) = s~^^'^, f2 = hi, /a = hj and f^is) = h{s)/s. Note that this is legitimate 
a p/n — pi > is sufficiently small; it yields the extra factor s^+^/p appearing below. 
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Upon using the pointwise bound f l6.2ip for the remaining k — 2 factors of components 
of /i.'s, (16.231) and the fact (since k <n) that 

we obtain the bound 

k 

|(F^,fc;/ii,...,/ifc+i)| <C|xrl^l J] (6.24) 

m=l 

By the induction hypothesis 

k 

m=l 

which together with (16.241) yields 

|(F^,,; h,..., Vi)| < C(a;)-l''l(a;)'=-('"-l^l)||/i||^, = C{xy-^\\h\\n.. (6.25) 

Suppose next that j = k + 1. Then we apply (I6.22p with m = 3, pi < p/n, p2 = p 
and p3 = q, fi{s) = s~"/^, /2 = hi and f3 = h yielding an extra factor s"/^. Upon using 
the pointwise bound (I6.2ip for the remaining k — 1 factors of components of /i.'s, (I6.23P 
and the fact that 

g-\ri\g\(:\g(k~m-l/p)gn/p <- 

we obtain again (I6.24p and whence (16.251) . 

Case Ail): 9''G(7(s)) = s-l'^^l(9^i9;^(,)G(sx + /t(s)); |Ci| = |C| - 1, H = k. Suppose 
first that j < k. Then we apply (16.221) with m = 3, pi < p/n, p2 = p and ps = q, 
/i(s) = /2 = hj and /3(s) = h(s)/s yielding an extra factor s^+^/p. Upon using 

(]6.2ip for the remaining k — 1 factors of components of /i.'s, (I6.23P and the fact that 

g-\v\g\(:i\g{k-l){l-l/p)gl+n/p 

we obtain again (I6.24p and whence (16.251) . 

Suppose next that j = k + 1. Then we apply (I6.22p with m = 2, pi < p/n, p2 = q, 
fi{s) = f2 = h yielding an extra factor s"'^^. Upon using (16.211) for the remaining 

k factors of components of /i.'s, (16.230 and the fact that 

g-lvlglCilgKT'-^/p) gn/p ^ 

we conclude as before. 

Case Aiii): d^Gi^is)) = s-K^\d^^^d;!^^^G{sx + K{s)y, IC2I = |C| - 2, \uj\ = k + 1. We 

apply (I6.22P with m = 2, pi < p/n, p2 = q, /i(s) = s""/^ and /2(s) = h{s)/s yielding 
an extra factor Upon using (I6.2ip for the remaining k factors of components of 

/i.'s, (I6.23P and the fact that 

g-\v\g\(2\gHl-l/p)^l+n/p 

we conclude as before. 

For the case B) we need (I6.23P with |?7| = |C| + ^, and we consider subcases Bi)-Bii) 
defined as before. 
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Case Bi): d''G{-f{s)) = s-\^\did';^^^-^G{sx + k{s)); \uj\ = k. Suppose first that j < k. 

Then we apply ( 16.22P with m = 3, pi < p/n, p2 = p and ps = q, = s""/^, /2 = hj 
and fsis) = h(s)/s yielding an extra factor gi+^/p. Upon using f l6.2ip for the remaining 
k — 1 factors of components of /i.'s, f l6.23p and the fact that 

we obtain the bound 

k 

\{F^,k;h,...,hk+i)\<C\x\'\x\-^^^ II \\hrr,\\\np\\h\\w. (6.26) 

m=l 

By the induction hypothesis 

k 

m=l 

which together with fl6.26p yields 

h,..., < C{xy-\^\{x)'-^--\<\^\\h\\u. = C{xy-^h\\n.. (6.27) 

Suppose next that j = k + 1. Then we apply ( I6.22p with m = 2, pi < p/n, p2 = q, 
/i(s) = s""/^, f2 = h yielding an extra factor s"/^. Upon using f l6.2ip for the remaining 
k factors of components of h,'s, fl6.23p and the fact that 

we conclude as before. 

Case Bii): 9''G(7(s)) = s-\<^\d^^'d^(^^^G{sx + /t(s)); |Ci| = |C| - 1, \io\ = k + 1. We 
apply fl6.22p with m = 2, pi < p/n, p2 = q, /i(s) = s""/^ and /2(s) = h{s)/s yielding 
an extra factor s^"*""/^. Upon using ( I6.2ip for the remaining k factors of components of 
h,'s, fl6.23p and the fact that 

g-\v\s\Ci\sH^-'i-/p)s'i-+n/p <- 

we conclude as before. 

For the case C) we need f l6.23p with |?7| = |^| + A; + 1. 

Case C): d'^G{-f{s)) = s-\'^\di&^^^^G{sx + k(s)); \uj\ = k + 1. 

We apply fl6.22p with m = 2, pi < p/n, p2 = q, fi{s) = and /2(s) = h{s)/s 

yielding an extra factor gi+^/p. Upon using fl6.2ip for the remaining k factors of com- 
ponents of /i.'s, fl6.23p and the fact that 

g-\v\g\C\gka-^/P)gl+n/p < 

we obtain the bound 

k 

\{F^,k;h,...,h+i)\<G\x\^\x\-\^\ n \\hj\\np\\h\\n^. (6.28) 

m=l 
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By the induction hypothesis 

k 



m=l 



which together with fl6.28p yields 

/ii, . . . , hk+i)\ < = C{xy-^h\\n.. (6.29) 

□ 

Remarks. 1) As noticed in the beginning of the above proof, using the proof of The- 
orem [T]l][ii]) we can improve (16.191) for |a| = 1 as 

\\d^4np < V^Cp for |a| = 1. (6.30) 

2) By integrating (I6.30p we obtain the bound 

\\i^\\np < V^Cp{x). (6.31) 

From the method of proof we have Cp — ?■ oo as p — ?■ oo. We remark that this feature 
is not an artifact of the proof. In fact for the example in Subsection 17. 31 the geodesies 
emanating from are rotating like logarithmic spirals showing that (I6.3ip is false for 
p = oo (here by definition = VFq^'°°(0, 1)'^). This partly explains why we worked 
with fl6.15bp rather than the e-independent pairing (16. 2p . 

Proof of Theorem \l.(A We mimic the proof of Proposition I3.1|[iv|) . Indeed combining 
(13. 7p and the bounds (16.190 and (I6.3ip (with p = 2) we show the following improvement 

of ^M) 

|<9"fi:(l)| < C(a;)i-I"l for all |a| < / - 1. 

Next we use again (13. 5p and obtain (I1.17ap . Obviously we can similarly use (13. 5p to 
prove ffLlTbl) . □ 



7. Examples 

In this section we present examples of metrics in the class O. Two of the examples 
are parameter-depending and are constructed by the exponential mapping from metrics 
that are not of order zero. 

7.1. Decaying potentials. Let be a radial function of class C' on R'^, l,d > 2, for 
which there are constants a > 0, A > 0, fi g]0, 2[ and a g]0, 2] such that 

-^(s)-^ < V{x) < -a{x)-^\ (7.1a) 

x-VV{x) + 2V{x) <aV{x), (7.1b) 

d'^Vix) = 0((x)-(^+l"l)) for |a| < /. (7.1c) 
Consider the functional J : M'^ x "H — )■ M given by 

J{x,k)= [' K{y{s)Ms)\'ds, (7.2) 
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where y{s) = sx + and K{x) = 2(A — V{x)) for A > 0, and consider the positive 
solution S{x) to the eikonal equation 

\VS{x)\^ = K{x) for X G \ {0}, (7.3) 

defined by 

Sir) = r y^Ki^dr, 
Jo 

or ahernatively by 

s\x) = m{j{x,k) -.Ken}. 

We introduce the diffeomorphism $ : R"' — )■ R'^ given by 

X 

^(x) = r(\x\)x = r(\x\)---, (7.4) 

fI 

where r is the inverse of the function 

p(r) = S(r)= / ^/K{T)dT. 



Note that ^{t^J K{{))x) = expo(tx), i.e. the exponential mapping at zero for the metric 
gij = K6ij of fl7.2p . A short calculation shows that if z{s) = ^{y{s)) and ijj{s) = 
yis)/\yis)\ then 

y-ujuj + —r-riy - y ■ 



--Py + —r-rP±y 



where P is the projection parallel to u and P± = I — P the projection onto {cu}"*". 
Therefore 

K{z)\z\' = \Py\' + f{\y\)\P^y\' 
= yG{y)y 

where 

G{y) = P + fi\y\)P^ (7.5a) 

and 

Thus, if y{s) = sx + k,{s) and z{s) = $(?/(s)) = s$(x) + h{s), with hi,h GTi, then 

J($(x),/i)= / y{s)G{yisms)ds=:E{x,K) 
Jo 

and therefore 

mf{E{x, K):KeH} = inf{ J($(x), h):heH} = S'^{^{x)). 
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Now we show that G satisfies (I1.3P with constants a and h independent of A > 0. It 
suffices to prove there exist c > and (7 > such that 




for all A > and r > 0. Note first that from (17. lap we have 

A + a(sr)-^ < ^ V{sr) ^ A + A{sr)-^' 
A + A(r)-/^ - \-V{r) - \ + a{r)-i' 

The lower bound follow from the fact that 



X + a{sr) ^ A + a(r) ^' a 



\ + A{r)-i^ ~ X + A{r)-^' ~ A 
for A > 0, r > and s G [0, 1]. Next we note that 



< ^ (7.6b) 



and therefore 




A + A(sr)-/^ , ^ 
ds < 




A + a(r)-'^ 



;1 + m2)m/4- 

The upper bound is obtained from the fact that the function 

is continuous for r > and satisfies 

lim0(r) = 1, and lim 0(r) = 2/(2 — yu). 

r—>0 r—^oo 

Thus, since it can easily be verified that G satisfies fll.2p . it follows that indeed G 
satisfies (11.30 . From from (I7.5ap it follows that G satisfies (I1.12p . A computation using 
fl7.1bp shows the bound (ll.lSp with any c < |/ sup /, whence G E O. 

Consequently, due to Theorem 11.41 if a real C'-potential on R'^, say V^, is sufficiently 
close to the radial potential V discussed above in the sense that for a sufficiently small 
e > 

\\V,-V\\i := supsup(x)l"l+^|9°(K(2;) -1/(|2;|))| < e, (7.7) 

\a\<l X 

then there exists a C'-solution to the eikonal equation (17. 3p with V V^. Indeed 
writing the solution from Theorem 11.41 defined in terms of the perturbed metric, say G^, 
as (in the coordinates y = $"^(2;)) we have S^{x) = 5e($~^(a;)) and we can use the 
estimates of Theorem 11.41 to compare the derivatives d°'Se of order |a| < 2 with the 
corresponding unperturbed quantities d°'S. The comparison in mind is given in terms 
of estimates exhibiting smallness in terms of the parameter e > 0. If we (for simplicity) 
assume that V is constant in a neighbourhood of zero then it is straightforward to show 
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that these estimates are uniform not only in > r for any r > but also in A > 0. 
All that is needed to show at this point is, cf. Remark II .51 121) . that 

sup IIG, - G\\i ^ as UK - V\\i 0. (7.8) 



A>0 



For (17. 8p we use in turn (and note for comparison) that the estimates of the constants 
in (II. 3p and (I1.13P given above are uniform in A > 0. The constructed solution has 
an application in scattering theory at low energies i.e. in the regime A — ?■ generalizing 
parts of |DS] . see 



7.2. Non-decaying potentials. Let V he a radial function of class C' on M'^, l,d > 2, 
for which there are constants a > 0, A > 0, /i g] — oo, 0] and a g]0, 2] such that 

< V{x) < -a(x)-'^, (7.9a) 

x-VV{x) + 2V{x) <aV{x), (7.9b) 

9"K(x) = 0((x)-(^+l"l)) for |a| < /. (7.9c) 

Notice that the these conditions are very similar to (I7.1ap -f l7.1cp . The only difference is 
that now /i < 0. We can again allow the inclusion of a positive energy, A > and obtain 
uniform estimates. In fact in this case (I7.9ap - fl7.9cp are invariant under the replacement 
V ^ V — X (up to a change of constants), so if not for the uniformity in A > we could 
in the following take A = 0. We get the bounds on / in a similar fashion. Note that we 
need to replace (I7.6ap and f l7.6bp by 

and 

< 4144:! < (7.10b) 

respectively. Again these estimates are for A > 0, r > and s e [0, 1]. We take the 
square root in (I7.10ap and ( ]7.10bp and integrate. As for (I7.10ap we then use a change 
of variables and obtain 




Obviously from (17.10bp we obtain 




V(r(fi)) 




We have verified that G given by (I7.5ap is a metric of order zero. As before we verify the 
bound (I1.13P with any c < f / sup /, whence G G C From this point we can proceed as 
before and introduce a class of perturbations K by (17. 7p (now with /i < 0) and indeed 
show the existence of a C'-solution to the eikonal equation (17. Sp with V ^ V^. 

We remark that the class of metrics discussed above by perturbing V{x) = —1/2 (in 
the particular case /i = A = 0) coincides with the class considered in |Ba] (here we also 
take / = 3). The parameter e > may play the role of inverse energy. In particular 
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the constructed solution to the eikonal equation was apphed in |ACH] in the study 
of scattering theory of Schrodinger operators with an order zero potential in the high 
energy regime. 

7.3. Order zero potential, logarithm orbits. We demonstrate in terms of an ex- 
ample from |HS] (see |HSt Example A. 6]) that perturbations of the Euclidean metric 
in the sense of this paper may involve somewhat exotic geodesies like logarithm orbits. 
This means that although the direction of any geodesic emanating from 0, say 7, and 
its velocity 7 are almost ligned up for all times, cf. Lemma 14.51 the geodesic is 
permanently rotating around as I7I — 00. 

Consider the symbol /i on x given hj h = h{x,C,) = g~^C,^, where the conformal 
(inverse) metric factor is specified in polar coordinates x = {rcos6,rsm6) as = 
g-se/fXj = f(^0 — elnr), X = x(r > 1). We assume / is a given smooth 27r-periodic 
function with max /' > 1 and that e > is small. Notice that the "x-space part" of the 
Hamiltonian orbits of this symbol are the geodesies of the metric Qij = gSij. Consider 
the orbits originating at (ro, 0; C, eC) where C > is arbitrary and tq > 2 is determined 
by the equation 

/'(^o) = (l + e')"'; ^^o = -elnro. (7.11) 
By assumption there exists at least one such solution. Let us assume that there are 
only a finite number of such solutions, say 6j, all being non-degenerate, f"{0j) 7^ 0. The 
x-space part of any corresponding orbit is the logarithmic spiral given by the equation 
9 — eliar = 9o- For f"{9o) > the orbit corresponds in reduced variables to a saddle, 
see |HS] . On the other hand for f"{9o) < the orbit corresponds to a sink. This means 
that generically the geodesies of the metric gij = g6ij emanating from are attracted 
to one of the logarithmic spirals associated to (17. lip and the condition f"{9o) < 0. 
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